Random sequential adsorption of partially oriented linear fc-mers on square lattice 
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Jamming phenomena on a square lattice are investigated for two different models of anisotropic 
random sequential adsorption (RSA) of linear /c-mers (particles occupying k adjacent adsorption sites 
along a line). The length of a k-mei varies from 2 to 128. Effect of k-mei alignment on the jamming 
threshold is examined. For completely ordered systems where all the /c-mers are aligned along 
one direction (e.g., vertical), the obtained simulation data are very close to the known analytical 
results for Id systems. In particular, the jamming threshold tends to the Renyi's Parking Constant 
for large k. In the other extreme case, when A;-mers are fully disordered, our results correspond 
to the published results for short /c-mers. It was observed that for partially oriented systems the 
jamming configurations consist of the blocks of vertically and horizontally oriented /c-mers (v- and 
/i-blocks, respectively) and large voids between them. The relative areas of different blocks and voids 
depend on the order parameter s, fc-mer length and type of the model. For small /c-mers (k < 4), 
denser configurations are observed in disordered systems as compared to those of completely ordered 
systems. However, longer fc-mers exhibit the opposite behavior. 

PACS numbers: 68.43.De, 68.35.Rh, 05.10.-a 

Keywords: jamming, Random Sequential Adsorption, Renyi's Parking Constant, critical exponent 



I. INTRODUCTION 

In the Random Sequential Adsorption (RSA), objects 
randomly deposit on a substrate; this process is irre- 
versible, and the newly placed objects cannot overlap or 
pass through the previously deposited ones |l| . The final 
state generated by the irreversible adsorption is a disor- 
dered one (known as the jamming state); no more objects 
can deposit in this state due to the absence of any free 
space of appropriate size and shape [l|-[3j . The fraction of 
the total surface occupied by the adsorbed objects, which 
is called the jamming concentration, is of central interest 
for understanding of the RSA processes. 

The RSA model studies of the objects with different 
shape (e.g., squares, ellipses, or stiff rods (needles)) have 
shown that jamming concentration depends strongly on 
the object shape and size (see, e.g., [1]). The RSA model 
has attracted a great interest as a tool for explanation 
of the wide class of irreversible phenomena observed in 
adsorption of chainlike and polymer molecules on ho- 
mogeneous and heterogeneous surfaces [4]. Jamming of 
short flexible linear chains on a square lattice has been 
studied by Becklehimer [5] and Wang [6]. Adsorption of 
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semi-flexible chains has been recently investigated using 
RSA model by Kondrat [7]. RSA of the binary mix- 
tures of extended objects (linear and bent chains) has 
been described by Loncarevic et al. [|[. RSA on a two- 
dimensional (2d) square lattice has been examined for the 
stiff and flexible polymer chains simulated by a sequence 
of lattice points forming needles, T-shaped objects and 
crosses, as well as flexible linear chains and star-branched 
chains consisting of three and four arms [9] . RSA model 
is widely used for simulation of thin film formation in the 
process of nanoparticle deposition [8|, llCHl3| . The inter- 
play between the jamming and percolation phenomena is 
discussed in several works (e.g., foL fl3l - [l5j ] ) . 

Important exact result for jamming concentration pj 
of /c-mers in one dimension (Id) random sequential ad- 
sorption has been reported by Krapivsky et al. [Hj] 



Pj = k / exp I —u — 2 

n \ .7 = 1 



1 - exp(-ju) 
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du. (1) 



2 « 0.86466 
*D(1) « 



In particular, the jamming concentration for dimers 
(k = 2) placed along one line is pj = 1 — e 
and for trimers (k = 3) is pj = 3D (2) — 3e 
0.82365, where D(x) = e~ x J* e £ dt is Dawson's inte- 
gral. For k — >• oo, the jamming threshold tends to Renyi's 
Parking Constant pj c R w 0.7475979202 (T3|. 

Many efforts have been concentrated on the study of 
deposition of /c-mers on a discrete 2d substrate. Most 
of previous works have been devoted to the deposition 
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of randomly oriented linear k-mers on square lattices. In 
square lattice models, the number of possible orientations 
of k-mers is restricted to 2. When orientation of k-mer is 
fixed, the situation is equivalent to the above mentioned 
case of deposition on a Id lattice. 

For lattice model, Manna and Svrakic observed that 
deposited k-mers (1 < k ^ 512) tend to align parallel to 
each other [Tsj|. They form large domains (stacks) with 
voids ranging from a single site up to the length of k-mer. 
The similar stacking and void formation is observed also 
for continuous RSA model of deposition of infinitely-thin 
line segments [19|, |20| . 

It is worth noting that principally different jamming 
behavior is observed for lattice and continuous models. 
For example, for deposition of extremely elongated ob- 
jects k — » oo (i.e., in the limit of infinite aspect ratio) 
the asymptotic jamming concentration Pj (oo) is in the 
off-lattice case [20| . The value of pj approaches zero with 
increase of the aspect ratio k and follows the power law 



Pj (k) = pj (oo) + a/k a 



(2) 



with Pj (oo) = and a = 1/(1 + 2y/2) « 0.26 0. 

In the discrete case, the asymptotic jamming concen- 
tration Pj (oo) remains finite [15, SiH, HI. In the latter 
case, the presence of finite coverage by the infinite k-mers 
has been inter pre ted as a consequence of the alignment 
constraint [2l|, [22|. 

Note that for deposition on a discrete lattice [HI, [H, 
l2ll [22I, the limiting jamming concentration Pj (00) de- 
pends upon the deposition mechanism. For the conven- 
tional RSA model (the vacant lattice site is randomly 
selected and any unsuccessful attempt of deposition of 
k-mer is rejected) and completely orient ationally disor- 
dered deposition of linear k-mers on square lattices, dif- 
ferent Monte Carlo studies have given the estimation 
Pj(oo) « 0.66 [H, SI. For, so called, "end-on" RSA 
model [2l| (in this model, once a vacant site has been 
found, the deposition is (randomly) attempted in all the 
directions until the segment is adsorbed or rejected), the 
Monte Carlo calculation has given the noticeably smaller 
value pj (00) = 0.583 ± 0.010 [Hj]. 

Different functions have been tested to fit the Pj(k) 
dependence for the deposition of randomly oriented k- 
mers on a two-dimensional square lattice (conventional 
RSA model) 0, [H, El EI]. MC data on jamming 
concentration Pj(k) obtained for line segments of length 
2 ^ k ^ 512 on square lattices of linear size L ^ 4096 
(preserving in all cases the ratio L/k > 8) have been 
fitted using the series expansion 



Pj(k) = Pj(oo) + a/k + b/k 2 



(3) 



Dj (OO) 



0.660 ± 0.002, a « 0.83 and 



and have giyenj 
b& -0.70 THlEl. 

The similar data for line segments of length 2 ^ k ^ 24 
on lattices of linear size L = 2000 have been approxi- 
mated using an empirical law [14j 



However, this equation is in fact a particular case of 
Eq. [3] under constrain of b = —a/2. It can be easily 
checked by using the substitutions of p* = Pj (00) + a/2 
and 7 = a I (2p*). Fitting of the numerical data presented 
in Table IJ of [14]) with Eq. Ogives the following estima- 
tions: Pj(oo) = 0.684±0.003, a = 0.59±0.01, p = 0.9998 
(coefficient of determination). 

Kondrat and Pekalski have reported the results for k- 
mers with length within 1 ^ k ^ 2000, lattice size L = 
30, 100, 300, 1000, 2500 and more than 100 independent 
runs [15]. Application of the power law (Eq. [2]) have 
given pj (oo) = 0.66 ± 0.01, a « 0.44 and b « 0.77. 

To our best knowledge, the very limited number of 
works have been devoted to jamming of non-randomly 
oriented linear k-mers and only particular case (k = 2) 
has been taken into consideration [23|, E3 . For isotropic 
problem, t he j amming concentration at k = 2 is pj ~ 
0.9068 0E|. 

In anisotropic problem, the vertical and horizontal ori- 
entations occur with different probabilities and degree of 
anisotropy can be characterized by the order parameter 
s defined as 



N\ + 



(5) 



p 3 (k)=p*(l- 1 (l-l/k) 2 ). 



(4) 



where N\ and N_ are the numbers of line segments ori- 
ented in the vertical and horizontal directions, respec- 
tively. 

Data of Monte Carlo simulations evidence that the 
value of pj decreases with order parameter s increase [23|, 
l2J]. In particular case of complete ordering, i.e., at s = 1, 
the problem becomes one-dimensional. An interesting 
finding is that in the limit of s — >• 1 the asymptotic frac- 
tion of dimers with horizontal direction does not vanish 
but equals to e~ 2 [l - exp(-2e- 2 )]/2 w 0.016046 @. 

The main goal of the present study is to investigate 
the effects of k-mei length, alignment, and the deposi- 
tion rules on the jamming threshold. This investigation 
is the natural development of the recent work devoted to 
the dimers, k = 2 [24]. This work discusses the jamming 
phenomena [l| for two different kinds of anisotropic se- 
quential deposition of linear k-meis on a square lattice 
(particles occupying k adjacent adsorption sites along a 
line). 

The paper is organized as follows. Section HI1 contains 
the basis of the two models of deposition of k-meis on a 
square lattice. The results obtained using finite size scal- 
ing theory and dependencies of the jamming threshold, 
p, and mean radius of pores, r, vs. order parameter, 8, 
are examined and discussed in details in Section HIB We 
discuss the dependence of the jamming threshold on the 
model parameters of interest in this Section, too. 
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Figure 1. Actual order parameter s, after each successful 
attempt to place a new k-mer (k = 3) for RSA and RRSA 
models. Lattice size is L = 729 and predetermined parameter 
is s = 0.8. 



II. DESCRIPTION OF MODELS AND DETAILS 
OF SIMULATIONS 

One can imagine suspension of linear line segments of 
length k (fc-mers) in the bulk volume under substrate. 
In our simulations, a square lattice of linear size L has 
been used as a substrate, and periodical boundary con- 
ditions have been applied both in horizontal and vertical 
directions (toroidal boundary conditions). The k-mers 
repulse each other and deposit one by one onto a sub- 
strate. The anisotropy of k-mer orientation in the sus- 
pension is predetermined and can be characterized by the 
input parameter s defined as in Eq. [5] 

The number of species in the suspension is supposed 
to be infinitely large, thus deposition doesn't change the 
anisotropy of the suspension. If different orientations of 
the deposited objects are not of equal probability (i.e., 
N\ 7^ N ), the definition of the jamming state is to be 
refined [24(. Let us assume N\ > N_. We define jamming 
for the fixed parameter s as a situation when there exists 
no possibility of deposition for any additional vertically 
oriented objects. Nevertheless, there may be places for 
horizontally oriented objects. 

Two different deposition models have been studied. 
The first is the conventional Random Sequential Adsorp- 
tion (RSA) model. In this model, the lattice site is ran- 
domly selected and an attempt of deposition of a k-mer 
with orientation defined by order parameter s is done. 
Any unsuccessful attempt is rejected and k-mer with a 
new orientation is selected. The second model, called by 
us Relaxation Random Sequential Adsorption (RRSA), 
is very similar to the RSA, however, the unsuccessful at- 
tempt is not rejected and a new lattice site is randomly 
selected until the object will be deposited. Note that in 
contrast to the known RSA model with diffusion (see, 
e.g., [26]), the RRSA model does not restrict the move- 



ment of species by the nearest sites only. The species may 
move all over the substrate searching for a sufficiently 
large empty space. In both models, the deposition ter- 
minates when a jamming state is reached along one of 
direction. 

Physically, the difference between RSA and RRSA 
models can reflect different binding of k-mers near the ad- 
sorbing substrate. In RSA model, binding of a k-mer by 
the substrate is weak and the k-mer returns to the bulk 
suspension after an unsuccessful attempt to precipitate. 
In RRSA model, binding of a k-mer to the adsorbing 
plane is strong, and k-mer has an additional possibility 
of joining the surface after an unsuccessful attempt. 

The differences between RSA and RRSA models can 
be evidently demonstrated by analysis of anisotropy of 
the deposited layer actually obtained in the course of 
adsorption (Fig. [1]). The preliminary study has shown 
that RSA model does not allow preservation of the or- 
der parameter s. In this model, the substrate "selects" 
the k-mer with appropriate orientation, and it results in 
deviation of predetermined order parameter s from the 
actually obtained one, so- The MC simulation evidences 
that the value of so noticeably exceeds the value of s. On 
the contrary, RRSA model better preserves the predeter- 
mined anisotropy, and sq « s (Fig. [TJ. 

In our study, the length of k-mers varies between 2 
and 128. To examine the finite size effect for RSA and 
RRSA models, different lattice sizes up to L = 100/c have 
been used. The number of runs varies from 10 to 1000 
depending of the lattice size L. 

Two different random number generators have been 
applied for filling in the lattice with objects (fc-mers) at 
given concentration and orientation. They are Mersenne 
Twister random number generator [28] with a period of 
2 19937 — 1 and the generator of Marsaglia et al. [29^ 
The results obtained using different generators are al- 
most undistinguishable. 

The connectivity of k-mers oriented in the vertical di- 
rection has been analyzed for jamming configurations, 
and the percolation threshold has been determined using 
the Hoshen-Kopelman algorithm (30| . 



III. RESULTS AND DISCUSSION 
A. Jamming configurations and its connectivity 

Jamming configurations obtained in the simulations 
are strongly dependent upon the order parameter and 
length of k-mers. However, visually, the difference in the 
structure of jamming patterns is not noticeable for the 
studied RSA and RRSA models. Examples of the typi- 
cal jamming configurations for RRSA model at different 
values of order parameter s are presented in Fig. [2j 

For randomly oriented linear k-mers (i.e., at s = 0) the 
typical domain structure in form of blocks of parallelly 
oriented k-mers has been observed. These blocks can be 



(a) RSA & RRSA, 5 = 0.0 



(b) RSA, s = 0.2 



(c) RSA, s = 0.6 




(d) RSA & RRSA, 5 = 1.0 (e) RRSA, s = 0.2 (f) RRSA, s = 0.6 



Figure 2. Examples of jamming configurations of /c-mers (k — 16) on a square lattice of size L = 2048 at different values of 
order parameter s. Vertical and horizontal orientations are represented by different grey levels in printed version and in red in 
blue in on-line version, empty sites are labeled white. Here, the part of the lattice sized 512 x 512 is shown. 



represented as the squares of size k x k [14|. One can 
present a jamming configuration as a combination of: 

• blocks of vertically oriented k-mers (^-blocks); 

• blocks of horizontally oriented k-meis (/i-blocks); 

• empty sites (voids). 

The relative area occupied by v- and ft-blocks is ap- 
proximately the same at s = (Fig. [2fa)). Increase of 
s results in decrease of the relative area occupied by h- 
blocks and ^-blocks become dominating in jamming pat- 
terns. Moreover, the relative area occupied by ^-blocks 
is visually larger for RRSA model than for RSA model 
(compare Fig. E(b),(c)) with Fig. E{e),(f)). Finally, at 
s = 1, the jamming configuration transfers into the ld- 
like domains (independent parallel Id jamming lines of 
k-meis) (Fig. [2fd) ) . 

It is interesting that the infinite connectivity between 
the similar v- or ^-blocks (i.e., percolation) is not ob- 
served for randomly oriented linear k-mers (at 8 = 0). 
The visual observations of jamming patterns show that 



connectivity between ^-blocks increases and between h- 
blocks decreases with increase of s. It can be easily 
demonstrated by analyzing the structure of the largest 
cluster of vertically oriented k-meis (Fig. [3]). At small s 
(e.g., at s = 0.05 in Fig.[3^a), (d)), the connectivity of v- 
blocks is limited and the largest cluster occupies only the 
finite part of the lattice. However, the size of the largest 
cluster is higher for RSA model than for RRSA model. 
Increase of s results in growth of the size of the largest 
cluster, and this cluster spans through the lattice at some 
threshold value of s c (see, e.g., Fig. [31(b) for RSA model 
and Fig. [31(f) for RRSA model). Examples of the perco- 
lation probability / vs. the order parameter s at k = 2, 8 
and L = 1024 are presented in Fig. 21 The threshold 
value of s c that corresponds to the percolation of verti- 
cally oriented k-meis is determined from the condition of 
/ = 0.5. 

The usual finite size scaling analysis of the percolation 
behavior is done, and it is obtained that s(L) follows 
scaling law governed by the universal scaling exponent v. 

\s c (L) - s c (oo)\ oc L~*, (6) 
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(d) RRSA, s = 0.050 (e) RRSA, s = 0.120 (f) RRSA, s = 0.245 

Figure 3. Examples of largest clusters of vertically oriented fc-mers (k — 16) for jamming configurations on a square lattice of 
size L — 2048 at different values of the order parameter s. 



where v = 4/3 is the critical exponent of correlation 
length for 2d random percolation problem [27j. 

For the thermodynamic limit (L — >• oo), dependencies 
of s c () vs. length of k-mei for RSA and RRSA models 
are presented in Fig. [5j the inset of this figure shows ex- 
amples of finite size scaling dependencies for k = 2 and 
k = 8. The observed percolation behavior of the ver- 
tically oriented /c-mers is rather different for RSA and 
RRSA model, which evidently reflects the difference in 
the structure of jamming patterns of these models. For 
the same /c, the RSA model gives lower value s c than 
RRSA model, e.g. k = 16, s c « 0.126 (RSA model) 
and s c w 0.240 (RRSA model). Note that in both mod- 
els for the dimers (k = 2), the threshold values of s c 
are rather close to « 0.21. However, the value of s c de- 
creases as the length of k-mei increases for RSA model, 
and opposite behavior is observed for RRSA model. Fi- 
nally, in continuous limit k — » oo, the difference between 
the threshold order parameters of RSA and RRSA mod- 
els becomes rather large, As c « 0.12. The difference in 
behavior of s c (k) observed for RSA and RRSA models 
evidences higher connectivity between v-blocks for RSA 
model than for RRSA model, and this tendency increases 



as the length of k-mer increases. 



B. Finite-size effects 

Commonly, it is assumed that finite-size effects on 
the jamming coverage of a periodic lattice decrease very 
rapidly with the lattice size L increase [3l|, [32[. For 
isotropic problem of adsorption of k-mers on a square 
lattice (5 = 0), it was observed that standard deviation 
of the coverage is a linear function oik/ L and it become 
negligible in the limit of k/L — >> [2lj . 

Our investigations show that in case of anisotropic 
problem the finite-size effects, as well as k and s, are 
rather sensitive to the type of adsorption model (RSA or 
RRSA). Fig. [6] demonstrates examples of jamming con- 
centration pj vs. s for both models and different lattice 
sizes L = 32-2048 in two particular cases when k — 3 and 
k = 8. For clearer demonstration of the finite-size effects 
on the jamming concentration, they are also represented 
in the form of dp = Pj(L — > 00) — Pj(L) vs. the inverse 
lattice size 1/L (see insets to Fig. [6]). 

In RSA model, (Fig. |B(a),(b) the finite size effects are 
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Table I. Estimated jamming concentrations, pj = Pj(L — >> oo), as a function of order parameter, s, at different values of k for 
RSA and RRSA models. Data Id correspond to exact results that are obtained by direct calculation using Eq|T] 
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Figure 4. Percolation probability / us. predetermined order 
parameter s for RSA and RRSA models at different values of 
k. Lattice size is L = 1024 and results are averaged over 100 
independent runs. 



rather small in the limits of s — >• and s —> 1. However, 
for the intermediate anisotropy s, they are large and dp = 
Pj(L — >• oo) — pj vs. 1/L looks like nonlinear even for 
the large lattice sizes L. From the other side, in RRSA 
model (Fig. |6fc),(d)) the finite size effects are very small 
for any k and s with the exception of the limit of isotropic 
problem, s ^ 0. 

The differences between the finite size effects of RSA 
and RRSA models at intermediate anisotropy (s = 0.8) 
and different values of k are presented in Fig. [7] in the 
form of pj vs. 1/L dependencies. For RSA model, these 
dependencies look like nonlinear even for large lattice 
sizes, and this fact obstructs application of the scaling 
relation in the thermodynamic limit (L —> oo) [2?) 

\Pj( L )-Pj(°°)\ ^L~i, (7) 

where v = 1.0 ±0.1 0- 

As a result, we have to utilize rather large lattices (L > 
512-2048) and nonlinear relations for p(l/L) in a form of 
polynomials of degree 3 for extrapolation of the results 
to thermodynamic limit (L —> oo) (see, e.g., Fig. [7]). For 
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Table II. Estimated parameters pj(oo), a and, a in a power law function Eq. [2] using for the fitting of a dependence of jamming 
concentration pj vs. length of linear segment k. The data are presented for RSA and RRSA models at different values of s. 
Data Id, exact, corresponds to fitting exact values that are obtained by direct calculation using Eq[T] In all cases the coefficient 
of determination p was higher than 0.9997. 
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Figure 5. Threshold order parameter s c = s c (L — » oo) vs. 
length of k-mer for RSA and RRSA models. Inset shows 
examples of scaling dependencies in the form of s c vs. L~ x ^ v 
for k — 2 and k — 8. Here, v — 4/3 is the critical exponent of 
the correlation length for 2d random percolation problem [23] • 



RSA model, this technique allows to get rather reliable 
estimations of jamming concentrations at L — >• oo with 
the error bar not exceeding ±0.002. 

For RRSA model, the finite size effects are very small 
at any L/k ^> 1 and s (with the exception of the limit of 
isotropic problem, s — »• 0); therefore, the final investiga- 
tions were performed for different k with L = 100/c and 
100 independent runs. Our estimation of the error bar 
for p(L —> oo) is about ±0.0001. For the particular case 
of the limit of isotropic problem, s — »• the procedure 
used for estimation of p(L — >• oo) is the same as for RSA 
model. 



Note that at moderate ratio L/k < 5 the finite 
size effects may be rather complex. The examples of 
Pj(L) — Pj (oo) vs. L/k dependencies for completely or- 
dered /c-mers, (s = 1) at different values of k are pre- 
sented in Fig. [8j). The observed noticeable oscillating 
scaling behavior evidently reflects commensurability of k 
and L values. 

The precision of the estimation of p(L — > oo) is tested 
for the limit cases of completely anisotropic (s = 1) and 
isotropic (s = 0) problems. 

The particular case of complete alignment (s = 1) cor- 
responds to the simplified Id problem when RSA and 
RRSA models are indistinguishable. The data presented 
in Table [I] show also that our results for both RSA and 
RRSA models at s = 1 are very close to the analyti- 
cal results calculated Id problem from Eq. (pQ). The nu- 
merically obtained data are precise within 4-5 significant 
digits. 

For isotropic problem (s = 0), both RSA and RRSA 
models give very close estimations of Pj(k) value. In this 
case, our results for pj for small k-mers (k = 2-8) (see 
Table |T]) are in a reasonable correspondence with the pre - 
viouslypublished data, e.g., 0.9068 p||, 0.9067(7) (pi, 
0.906 [14] (k = 2), 0.8465 0, 0.847 pj|" (k = 3), 0.811 Hi 
(k = 4) and 0.757 [14] (k = 8). 

The similar jamming behavior of the RSA and RRSA 
models in the isotropic case s = is expected. For RSA 
model in the isotropic case, rejection of the unsuccess- 
ful attempt is followed by the next choice of k-mers with 
random orientation. In this case, the RSA model also 
allows preservation of the predetermined order parame- 
ter. That is a reason why the RSA and RRSA models 
give the similar estimations for the values of Pj(L —> oo). 
However, the amplitude of order parameter fluctuations 
during the deposition of k-mers may be rather different 
for RSA and RRSA models, and in fact it results in a 
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Figure 6. Jamming concentration p vs. order parameter s at different values of L. Insets show dp = Pj(L — >> oo) — pj vs. 1/L 
for different values of s. 



noticeable difference of the finite-size scaling effects at 
5 = observed for RSA and RRSA models (see Fig. [6j). 

Finally, the data on p(L — » oo) vs. s dependencies for 
RSA and RSSA models at different k obtained in a result 
of the described scaling analysis are presented in Table HI 



C. Jamming concentration 

1. Jamming concentration vs. order parameter 

Fig. [9] demonstrates behavior of jamming concentra- 
tions, pj = Pj(L —> oo) , as a function of order parameter, 
5, for RSA and RRSA models, respectively. 

The Pj(s) dependencies are rather different for short 
and long /c-mers. For instance, if k-mers are short 
(k ^ 4), more dense configurations are observed for disor- 



dered systems (s = 0) as compared to those of completely 
ordered case (s = 1). However, an opposite behavior is 
observed for longer /c-mers. For RSA model, the value 
of Pj(k) monotonically decreases (k ^ 4) or increases 
(k > 4) when the order parameter s increases. For RRSA 
model, the value of Pj(k) goes through the minimum at 
k ^ 3 (RRSA-model) when s increases. 

At each value of /c, the jamming concentrations pj are 
practically the same at s = and s = 1 for the both RSA 
and RRSA models. However, in the intermediate region 
of s we observe £>j(s)(RSA)> pj(s)(RRSA). The maxi- 
mal difference between pj(s) (RSA) and £>j(s)(RRSA) is 
observed at s « 0.5-0.7. 

Such behavior may be explained as follows. In RSA 
model, the substrate "selects" the k-mer with appropri- 
ate orientation. It results in amplification of the actual 
order parameter so > s (Fig. [TJ. 
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Figure 7. Jamming concentration pj vs. 1/L for different values of k = 3(a), 8(b), 64(c) and fixed order parameter, s = 0.8. 
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Figure 8. Difference of jamming concentrations Pj(L) —pj(oo) 
vs. L/k ratio for completely ordered /c-mers, (s = 1) at dif- 
ferent values of fc(RSA and RRSA models). 



Such amplification reflects intensive nucleation of v- 
domains at initial stages of deposition. At s > 0, the 
^-domains merge into percolating clusters when s > s c . 
For well oriented systems (at s ~ 0.5 — 0.7), the jamming 
structures in RSA model are large ld-like ^-domains with 
small inclusions of /i-domains and voids. 

In RRSA model, the order parameter keeps the level 
of approximately so ~ s (Fig. [T]). The packing of fc-mers 
in the form of fo-and ^-domains arising at initial stages 
of RRSA deposition is sparser than for RSA model. As a 
result, the ^-domains merge into percolating clusters at 
higher level of s c in RRSA model than in RSA model. 
Finally, for the well oriented systems (at s « 0.5-0.7), 
the jamming structures of RRSA model include higher 
number of /i-domains with large voids between them than 
those of RSA model. 



The effective pore radius r in the jammed structures 
is calculated as 



(8) 



where 1 — p j is the specific area of pores, a = A/L 2 is 
the specific total perimeter of pores, and A is the total 
perimeter of pores. 

The same as for jamming concentration p^ the values 
of r are practically equivalent at s = and s = 1 for both 
RSA and RRSA models. Fig. [10] compares the effective 
pore radius r vs. value of k (a) and jamming concentra- 
tion pj = Pj(L — > oo ) (b) at s = and 8 = 1 (RSA and 
RRSA models). The value of r increases with k and s 
values growth (Fig. [Tola)). Moreover, direct correlations 
between the pore radius r and the jamming concentration 
Pj are observed. Increase of the jamming concentration 
evidently reflects decrease of the effective pore radius r 
(Fig. Mb)). 

Interesting correlations are also observed between 
the effective pore radius r and the order parameter s 
(Fig. ITT]). In close analogy with behavior of the jamming 
concentration pj( s ), dependencies r were rather different 
for short (k < 4) and longer /c-mers. If /c-mers are short, 
increase of s results in increase of the effective pore radius 
r. The opposite behavior is observed for longer /c-mers 
(k > 8). At same values of s and fc, the RSA model 
demonstrates better jamming packing and smaller effec- 
tive pore radius. 



2. Jamming concentration vs. length of k-mer 

For Id problem, the jamming concentration vs. length 
of k-mer can be evaluated exactly as expressed by Eq. [TJ 
Surprisingly enough, this very complex Pj(k) dependence 
may be fitted rather well by a simple power law (Eq. [2]) 
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Figure 10. Effective pore radius r vs. value of k (a) and jamming concentration pj = Pj(L — » oo) (b) at s = and s = 1 (RSA 
and RRSA models). 



with ft (oo) = c R [17], a = 0.2277 ± 0.0019, a = 1.0111 ± 
0.0023, p = 0.99997 (coefficient of determination). 

Fig. [T2l shows jamming concentration us. length of 
linear segment k at different anisotropy s for RSA and 
RRSA models). It has been found that numerical results 
may be rather well fitted by a power law function Eq. [2j 
However, for problem of nonrandomly oriented linear k- 
mers, pj(oo), a and, a depend on the order parameter 
s (Table HI]). For instance, for isotropic problem (s = 0) 
the obtained data (Table HI]) are in satisfactory agreement 
with the previously reported data: Pj(oo) = 0.66 ± 0.01, 
a ~ 0.44 and a w 0.77 [15|- The fitting parameters ob- 
tained for completely ordered /c-mers, (s = 1) (Table HI]) 



are also in satisfactory correspondence with the theoret- 
ical predictions for Id problem. 

We should emphasize that for the studied problem of 
nonrandomly oriented /c-mers, the jamming concentra- 
tion decreases as an inverse power of the linear segment 
length, i.e., pj ex k~ a and the exponent a is not univer- 
sal. Note that power law behavior with exponent a ~ 0.2 
was observed for off-lattice RSA of rectangles [19| and it 
was suggested that indicates about fractal structure of 
the jamming networks with a dimension df = 2 — a ~ 0.2. 
For completely ordered /c-mers, (s = 1) the value of df is 
close to 1 and it is expected for Id problem. For disor- 
dered systems at s < 1 the fractal dimension df increases 
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Figure 11. Relative effective pore radius r/ro vs. order pa- 
rameter s at different values of k for RSA and RRSA models). 
Here, r Q — r(s — 0). 



Figure 13. Fractal dimension of the RSA and RRSA jamming 
networks df vs. order parameter s. 
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Figure 12. Jamming concentration^ vs. k at different values 
of order parameter s. Inset shows dp = pj — Pj(k — > oo) 
vs. k, the different slopes correspond to the different scaling 
exponents a. (RSA and RRSA models). 



for RRSA model and goes through minimum for RSA 
model when the order parameter s decreases (Fig. [T3|) . 



IV. CONCLUSION 

Two different models describing jamming of the par- 
tially ordered linear k-mers (k = 2-128) with predeter- 
mined order parameter s on a square lattice are discussed. 
In usual RSA model, the substrate "selects" the k-mer 
with appropriate orientation, which results in amplifi- 
cation of the actual order parameter, so > s. In re- 
laxation RSA model (RRSA), the order parameter re- 
mains at the predetermined level, so ~ s. The similar 
jamming behavior is observed at different values of k 



for both RSA and RRSA models in problems of com- 
pletely disordered (at s = 0) and ordered (at s = 1) 
systems. Our new simulation results for jamming con- 
centrations pj are in excellent agreement with simula- 
tion data (s = 0, k = 2,3,4,8) [fl, [3, GJ, 0, [H and 
analytical results [16| previously published for the prob- 
lem of completely disordered system. For short k-meis 
(k ^ 4), more dense configurations are observed in disor- 
dered systems (s = 0) as compared with those completely 
ordered (s = 1). However, an opposite behavior is ob- 
served for longer fc-mers. In partially oriented systems 
(s < 1), the jamming configurations consist of the blocks 
of vertically and horizontally oriented k-meis (v-and h- 
blocks, respectively), and large voids between them. The 
^-blocks merge into the percolation cluster when the or- 
der parameter s exceeds some critical value s c , which 
depends upon k and model of deposition. The critical 
exponent of random percolation v = 4/3 is observed in 
the scaling relations. It is demonstrated that in the in- 
termediate region of < s < 1 the RSA process allows 
obtaining of better jamming packing and smaller effec- 
tive pore radius than RRSA process at the same values 
of s and k. For RSA model, the value of Pj(k) monoton- 
ically decreases (k ^ 4) or increases (k > 4) when the 
order parameter s increases. For RRSA model, the value 
of pj (k) goes through the minimum at k ^ 3 when s in- 
creases. Finally, it is found that numerical results for pj 
vs. k dependence may be rather well fitted by a simple 
power law function pj ex k~ a with parameters dependent 
on the order parameter s. The power exponent a grows 
when the order parameter s increases from a w 0.72 at 
5 = to a « 1 at s = 1. The jamming networks display 
fractal properties with fractal dimension df w 1 for com- 
pletely ordered /c-mers, (s = 1) and value of df increases 
for RRSA model and goes through minimum for RSA 
model when the order parameter s decreases (Fig. [T3|) . 
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